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We analytically show that it is possible to perform coherent imaging by using the classical cor- 
relation of two beams obtained by splitting incoherent thermal radiation. A formal analogy is 
demonstrated between two such classically correlated beams and two entangled beams produced by 
parametric down-conversion. Because of this analogy, the classical beams can mimic qualitatively 
all the imaging properties of the entangled beams, even in ways which up to now were not believed 
possible. A key feature is that these classical beams are spatially correlated both in the near-field 
and in the far-field. Using realistic numerical simulations the performances of a quasi-thermal and a 
parametric down-conversion source are shown to be closely similar, both for what concerns the reso- 
lution and statistical properties. The results of this paper provide a new scenario for the discussion 
fN| ' of what role the entanglement plays in correlated imaging. 

>v 

PACS numbers: 42.50.Dv, 42.50-p,42.50.Ar 

I. INTRODUCTION 

The topic of entangled imaging has attracted noteworthy attention in recent years 0,|2)|3,0)IEIEI3>I3)I3- This 
technique exploits the quantum entanglement of the state generated by parametric down-conversion (PDC), in order 
to retrieve information about an unknown object. In the regime of single photon-pair production of PDC, the photons 
of a pair are spatially separated and each propagates through a distinct imaging system, usually called the test and 
the reference arm. An object is located in the test arm. Information about the spatial distribution of the object is 
I/"") | obtained by registering the coincidence counts as a function of the transverse position of the photon in the reference 
arm, which holds a known reference system 0, 0, 0, 0, 0, Ej ■ In the regime of a large number of photon pairs, this 
procedure is generalized to the measurement of the signal-idler spatial correlation function of intensity fluctuations 
0- Such a two-arm configuration provides more flexibility in comparison with standard imaging procedures. For 
example, there is the possibility of illuminating the object at a given light frequency in the test arm and of performing 
Qh! a spatially resolved detection in the other arm with a different light frequency, or of processing the information from 
+1> ■ the object by only operating on the imaging system of the reference arm 6]. In addition, it opens the possibility for 
performing coherent imaging by using, in a sense, spatially incoherent light, since each of the two down-converted 
beams taken separately is described by a thermal- like mixture and only the two-beam state is pure (see and 0). 

In this paper we show that such a scheme can be implemented using truly incoherent light, as the radiation produced 
by a thermal (or quasi-thermal) source. A comparison between thermal and photon-pair emission was performed in 
[l0(, where an underlying duality accompanies the mathematical similarity between the two cases. Here we consider 
a different scheme (Fig. ^) appropriate for correlated imaging, in which a thermal beam is divided by a beam-splitter 
(BS) and the two outgoing beams are handled in the same way as the PDC beams in entangled imaging. Our analysis 
points out a precise formal analogy between the PDC and the thermal case. This analogy opens the possibility 
for using classically correlated thermal light for correlated imaging in the same way as entangled beams from PDC. 
Currentlythere is a very lively debate whether quantum entanglement is a necessary ingredient to perform correlated 
imaging [E 11 • The first coincidence imaging experiments were performed using entangled photons from PDC 
0,3- At that time, the Authors of Ref. Q suggested that "it is possible to imagine some type of classical source that 
could partially emulate this behavior". A more recent theoretical analysis gave arguments that "the distributed 
quantum-imaging scheme truly requires entanglement in the source and cannot be achieved using a classical source 
with correlations but without entanglement" . The topic became hot after the ghost image experiment of Ref. was 
successfully reproduced using classically correlated beams p|. In this experiment a classical source produced pairs 
of single-mode angularly correlated pulses that served as classical analogs of momentum correlated pairs of photons 
produced by PDC. In the accompanying theoretical discussion, the Authors presented arguments that although the 
results of any single experiment in quantum imaging could be reproduced by classical sources with proper statistical 
correlation, a given classical source cannot emulate the behaviour of a quantum entangled source for any arbitrary test 
and reference systems. In agreement with this we showed in a recent theoretical Letter that on the one hand the 
results of any single quantum-imaging experiment performed with entangled beams in a pure state could be exactly 
reproduced by using separable (i.e. non-entangled) states. On the other hand, we showed that a key feature of the 
entangled state produced by PDC is the simultaneous presence of spatial correlation at the quantum level in the 
near field and the far field of the signal- idler beams (see also [Tl|V This corresponds to the simultaneous quantum 
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FIG. 1: Correlated imaging with incoherent thermal light. The thermal beam a at the beam-splitter BS is divided into two 
beams, 61 and 62, which travel through respectively a test and a reference system, described by their impulse response functions 
hi and hi- The test arm 1 includes an object. Detector D\ is either a point-like detector or a bucket detector. D2 is an array 
of pixel detectors, v is a vacuum field. 



correlation of position and momentum of the photons in each pair. We showed that this feature could be exploited to 
produce both the image and the diffraction pattern of an object by solely operating on the reference arm, and argued 
that this could be done only in the presence of quantum entanglement. This interpretation was received rather well 
in the quantum imaging community and was generally viewed as a possibility to discriminate between the presence 
of quantum entanglement and classical correlation in the source. 

In this paper we will analyse a counter-example, which partially contradicts the picture emerging from Refs. 
The classical, thermal beams created by splitting thermal radiation have several features that distinguish them from 
other non-entangled beams considered in the previous literature 0,13- First of all, they are spatially multi-mode, 
unlike those considered in Ref. Second, as we will see in Section IIVI they have an imperfect level of spatial 

correlation, limited by the classical shot-noise introduced by the beam-splitter. Nonetheless, they exhibit a high level 
of spatial correlation in both the near-field and far-field planes, in contrast to the separable mixtures considered in 
Ref. 7] that were perfectly correlated in either plane. We will show that this imperfect spatial correlation in both 
planes, although being completely classical, is enough to qualitatively reproduce all the features of the entangled 
imaging, provided that the spatial coherence properties of the thermal source are properly engineered. Finally, they 
are probably the best classical analogs of the entangled signal-idler beams produced by PDC, since the marginal 
statistics of the signal or idler beam alone is a thermal statistics. Thus, they should provide an optimal test bed for 
understanding the role of entanglement in correlated imaging. 

In Section II we demonstrate theoretically the analogy between thermal and entangled PDC beams in correlated 
imaging. Section III discusses in analytical terms a specific imaging scheme. In Sec. IV we discuss the origin of the 
spatial correlation in the thermal case and relate it to the entangled case. In Sec.|3the performances of the two cases 
are compared and we show a key numerical example that drops the spatio-temporal translational invariance assumed 
in the analytical treatment. Section VI contains the conclusions and a discussion. Finally, the Appendix discusses the 
degree and visibility of the correlation between the thermal beams when the finite detection area of the measurement 
apparatus is taken into account. 



II. ANALOGY BETWEEN THERMAL AND ENTANGLED BEAMS IN CORRELATED IMAGING 

In this section we are going to show a close analogy between the use of thermal light and entangled beams from 
PDC in the imaging schemes based on correlation measurements. For the sake of comparison, the two cases will 
be treated in parallel. In the analytical treatment we consider for simplicity only spatial variables and ignore the 
time argument, which corresponds to using a narrow frequency filter. We will come back to this point in Sec.[V] In 
addition, we assume translational invariance in the transverse plane, which amounts to requiring that the cross-section 
of the source is much larger than the object and all the optical elements. 

In the entangled case, the signal and idler fields are generated in a type II crystal by a PDC process. Our 
starting point is the input-output relations of the crystal, which in the plane- wave pump approximation read [tI IIIIH^ ] 



bi(q) = Uiidoitf) + Vi{q)a]{-q) i^j = 1,2 



(1) 
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Here, bi(q) = / |f e~ 1 ^' x b i {x), where bi(x) are the signal (i — 1) and idler (i = 2) field envelope operators at the 
output face of the crystal (distinguished by their orthogonal polarizations), x being the position in the transverse 
plane. <Zj , i = 1, 2 are the corresponding fields at the input face of the crystal, and are taken to be in the vacuum 
state. The gain functions Ui, Vi are for example given in 

In the thermal case, we start from the input-output relations of a beam splitter 

bi(x) — ra(x) + tv(x) , 62 (x) = ta(x) + rv{x) , (2) 

where t and r are the complex transmission and reflection coefficients of the mirror, a is a thermal field and v is 
a vacuum field uncorrelated from a. We assume that the thermal state a(x) is characterized by a Gaussian field 
statistics, in which any correlation function of arbitrary order is expressed via the second order correlation function 
14]: 

T(x,x') = (a + (x)a(x')> 

^ e -»*<*-*')<„($)>th. (3) 

Here {n(q))th denotes the expectation value of the photon number in mode q in the thermal state. In writing the 
second line of this equation, we implicitly used the hypothesis of translational invariance of the source, under which 
T(x, x') = T(x — x'). In particular, the following factorization property holds [l4|: 

(: a\x)a{x') a\x")a{x"') :) = (a^x)a(x'))(a^ {x")a(x"')) + { a\x)a(x' ")){a\x ")a{x") ) , (4) 

where : : indicates normal ordering. 

In both the PDC and the thermal case each of the two outgoing beams travels through a distinct imaging system, 
described by its impulse response functions h\{x\,x[) and ^-2(^2, ^2); respectively (see Fig. 0. Arm 1 includes an 
object. Beam 1 is detected by D\, which is either a point-like detector or by a "bucket" detector which collects all the 
light in the detection plane [g; in any case D\ gives no information on the object spatial distribution. In the other 
arm the detector D2 spatially resolves the light fluctuations, as for example an array of pixel detectors. The fields at 
the detection planes are given by 

Ci{xi) = dx-hi(xi,Xi')bi(xi') + Li(xi) £ = 1,2, (5) 



where £1,^2 account for possible losses in the imaging systems, and depend on vacuum field operators uncorrelated 
from b\, 62 ■ Information about the object is extracted by measuring the spatial correlation function of the intensities 
detected by D\ and D2, as a function of the position X2 of the pixel of D2: 

{h(xi)h{x2)) = {c\{xi)ci{xi)c\{x2)c 2 {x2)) ■ (6) 
All the object information is concentrated in the correlation function of intensity fluctuations: 

G(x u x 2 ) = {h{Sx)h{x2)) - {h{Sx)){Hx2)) , (7) 

where (Ii(x*i)) = (c\[xi)ci(xi)) is the mean intensity of the i-th beam. When using a bucket detector in arm 1, the 
measured quantity corresponds to the integral over X\ of both sides of Eq. Q. Since c\ and c 2 commute, all the terms 
in Eqs. ©,0 are normally ordered and Li,L 2 can be neglected, thus obtaining 

G(xi,X2)= J dx[ J dx" J dx 2 J dx^h^xi, x")h\(x\, x[)h2(x2, xl 1 r )h2(x 2l x'^) 

(blixi'Mx^bUxMxi)) - (blixi'Mxi)) (btiZ2)b2&))] . (8) 

In the thermal case, by taking into account the transformation @ and that v is in the vacuum state, 61 and &2 in 
Eq. © can be simply replaced by ra and ta, respectively. Next, by using Eq. we arrive at the final result 



G t h(^i,^2) = \trf 



2 

dx[ I dx2h\(xi,x[)h 2 {x2,X2)(a^ {x[)a(x2j) 



(9) 



where (a^ (x[)a(x2)) 1S given by iJSJl. 
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Similarly, the four-point correlation function in Eq. 
As it can be obtained from Eq. (see also 



has special factorization properties also in the PDC case. 



{b{{^)b x {x[)b\{x^{xi)) = {b\^)b x {x[)){b\{S!l)b 2 ^)) 

+ (bl(Zl')bU^'))(bi(Z{)b2&)) 

By inserting this result in Eq. one obtains 



G P nc(xi,x 2 ) 



dx{ / dx 2 hi(xi, x[)h 2 (x 2 , x 2 )(bi (x[)b 2 (x 2 )) 



where by using relations 



(h(x{)b 2 (x^)) 



^- S ^Ui(q)V 2 (-q) 



(10) 



(11) 



(12) 



At this point the analogy between the results in the two cases clearly emerges. Apart from the numerical factor 
\tr\ 2 and the presence of hi instead of hi, the thermal second-order correlation (aJ (x)a(x')) in Eq. 10 plays the same 
role as the PDC signal-idler correlation function (bi(x)b 2 (x')) in Eq. Ijlip. Consequently from Eqs. © and 1|12|) . the 
thermal mean photon number {n(q)) t h plays the same role as Ui(q)V 2 {— q) in the PDC case. The correlation function 
{a^ (x)a(x')) governs the properties of spatial coherence of the thermal source The correlation length, or 

transverse coherence length Z co h, is determined by the inverse of the bandwidth Aq of the function (n(q))th- The 
same comments hold for the correlation (bi(x)b 2 (x')), and the function Ui{q\V2(— q) in the entangled case. Most 
importantly, unlike the results for the separable states considered in Refs. yj Q, in both Eqs. © and Ijllfl the 
modulus is outside the integral; this is the feature that ensures the possibility of coherent imaging via the correlation 
function (see, e.g., 



III. IMAGING SCHEMES: IMAGE AND DIFFRACTION PATTERN OF AN OBJECT 

Let us now analyse two paradigmatic examples of imaging systems, borrowed from the discussion of and sketched 
in Fig. [3 In both examples the setup of arm 1 is fixed, and consists of an object, described by a complex transmission 




FIG. 2: Imaging scheme. L denotes two identical lenses of focal length /. Di is a point-like detector. The distance z is either 
z = / or z = 2/. 

function T(x), and a lens located at a focal distance / from the object and from the detection plane. Hence, 

hi(xi,x[) = ---j cxp f-^-Xi ■ x{\ T{x[) , (13) 

with A being the wavelength. In arm 2 there is a single lens placed at a distance z both from the source and from the 
detection plane 2; for simplicity we take the two lenses identical. 
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In the first example we assume z = f so that h 2 (x 2 , x 2 ) = — TJ ex P ("Tf^ 2 ' ^2) ' ^ msei "ting these propagators 
into Eq. (O and taking into account Eq. |J3}, we obtain 



G t h(^l,^2) oc 



. . _ 2tt.. ~ _ . 2vr\ 

(n(-a!2^-)}th-r I {x 2 ~ x Vjl J 



(14) 



where T(q) = J ^ e -^ s T( x) is the amplitude of the diffraction pattern from the object. This has to be compared 
with the result of the entangled case (see Eq. (7) of 0), 

Gpoc(xi,X2) oc 



Ux{S 2 — )V 2 (X 2J J )T [{X2 + J 



(15) 



where the combination x 2 + x\ appears instead of X2 — Xi and J7iV2 instead of (n)th- In both the thermal and 
the PDC case the whole diffraction pattern from the object can be reconstructed via the correlation function. This 
holds provided that the spatial bandwidth Aq is larger than the maximal transverse wave-number q in the diffraction 
pattern, or equivalently, provided that Z C oh < Iq, where l is the smallest scale of variation of the object spatial 
distribution. Thus both cases have best performances of the scheme when spatially incoherent light (l co h ~ * 0) is 
used. In contrast, as it is well known, when ^ co h < l no information about the diffraction pattern of the object can 
be obtained without the correlations, i.e. if we detect the light intensity distribution in arm 1 with an array of pixels. 
In fact, one can easily obtain that 



dq 



f(x 



2tt 

aT 



9) 



(n(q)) 



th 



(16) 



For Z co h < Z , (n(q))th can be taken out of the integral, and the resulting expression does not depend on x\ any more. 
We incidentally remark that the result of Eq. (|14|l differs from what one would obtain with a standard Hambury-Brown 
and Twiss scheme [T^ |. where the object is placed in the thermal beam a before the beam-splitter. In that case, one 
would retrieve the Fourier transform of the modulus square of the object transmission function (see e.g. p^l. thus 
losing any phase information. In our scheme instead, where the object is located in only one arm of the two, phase 
information about the object can be extracted and e.g. the diffraction pattern from a pure phase object can be 
reconstructed. 



In the second example, we set z — 2f, so that h,2(x2,x 
and taking into account (|13|) . we get: 



= S(x 2 + x'2) ex P ( ~ i|^2| 2 ^j ) ■ Inserting this in Eq. © 



Gt-h{xi,x 2 ) cx 



dx[T(x[ + x 2 )T* {x[)e^- 



\T(-x 2 ) 



(17) 



(18) 



where in the second line Z co h < l Q was assumed. Since the correlation function T{x[ J rX 2 ) — (<v(x{)a(— X2)) is non-zero 
in a region of size Z co h around x{ — —x 2 , this condition ensures that T(x[) is roughly constant in this region and it 
can be taken out from the convolution integral in |17jl. thus obtaining (|18|l . In this example the intensity correlation 
function provides information about the image of the object. In the general case <|17[) . the image reconstructed via 
the correlation function is a convolution of the object image with the second order correlation function J3J); therefore 
the thermal coherence length Z co h fixes the resolution of the imaging scheme. 

Under the same assumption Z co h < l Q , a similar result holds for the case of entangled beams (see Eq. (8) of 0)[25| 
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{xi,x 2 ) oc \T(-x 2 )\ 
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Vo - 
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(19) 



Also in this case the resolution of the scheme is limited by the finite transverse coherence length of the PDC beams. 

This section shows that the classical correlation of the thermal beams offers imaging capabilities similar to those 
of the entangled PDC beams; both the image and the diffraction pattern of an object can be reconstructed and we 
can pass from one to the other by only operating on the optical setup in the reference arm. The performances of 
the imaging schemes based on classical correlation and on entanglement will be compared in Sec. both for what 
concerns the visibility, the statistical properties as well as the spatial resolution. The next section, instead, will explain 
the basic mechanism that allows correlated thermal imaging with such a degree of flexibility. 
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IV. NEAR AND FAR-FIELD CORRELATION IN THE THERMAL AND ENTANGLED CASE 

As explained in detail in Ref. |jj the imaging schemes described in the previous section have a peculiar feature. In 
the z = f scheme the diffraction pattern reconstruction is made possible by the presence of spatial correlation in the 
far field of the correlated beams (momentum correlation of the photons). In the z — 2/ scheme, it is the presence of 
spatial correlation in the near-field (position correlation of the photons) that ensures the possibility of reconstructing 
the image. Our results for the thermal case may hence appear surprising if one has in mind the case of a coherent 
beam impinging on a beam splitter, where the two outgoing fields are uncorrelated, i.e. G{S\,X2) = 0. However, 
when the input field is an intense thermal beam, i.e. the photon number per mode is not too small, the two outgoing 
beams are well correlated in space both in the near-field and in the far-field planes. 

To prove this point, let us consider the number of photons detected in two small identical portions S ("pixels") 
of the thermal beams in the near field immediately after the beam splitter, Ni = J s dxb\(x)bi(x) i = 1,2, and the 
difference iV_ = Ni — N 2 . Making use of the transformation JS}, for |r| 2 = |i| 2 = 1/2 it can be proven that the 
variance (SN^_) = (N^_) — (-/V_) 2 is given by 

(SNl) = (Ni) + (JV a ) , (20) 

which corresponds exactly to the shot noise level. Remarkably Eq. Q20(l holds regardless of the statistical properties 
of the input beam a provided that in the other input port there is the vacuum. On the other hand, by using the 
identity (SNl } = (8N?) + (SNg) - 2(6N 1 5N 2 ) and taking into account that (5N?) = {8N%) f or | ? f = ^ the degree 
of spatial correlation is described by 

c def (SNxSN^ = _ (Nj) 

For any state < \C\ < 1, where the upper bound is imposed by the Cauchy-Schwarz inequality. The lower bound 
corresponds to the coherent state, for which (SNf) — (Ni) . For the thermal state, there is always some excess noise 
with respect to the coherent state (SNf) > (Ni), so that the correlation 121|) never vanishes. Remarkably, a high 
degree of spatial correlation between beams b\ and b 2 is ensured by the presence of a high level of excess noise in the 
input beam. As shown in detail in the Appendix, for thermal systems with a large number of photons, provided that 
the pixel size is on the order of Z co h or larger, (Ni)/(SNi) <C 1, and C can be made close to its maximum value. 

Even more important, in the absence of losses it is not difficult to show that Eqs. I|2U|) and l|21|l hold in any plane 
linked to the near- field plane by a Fresnel transformation. Let us assume that the propagation of beams 61,62 is 
described by a linear and unitary kernel H, bn,i{x) — J dx'H(x, x')bi(x'), i = 1,2 . Then the form of the beam- 
splitter transformation (|2J) is preserved during propagation, provided that the thermal field a is substituted by the 
propagated field a#(x) = J dx'H(x, x')a(x'). Hence Eqs. <j2*U|) and (|2*T)> also hold for bn,i and bu,2t because these 
equations are just a consequence of the beam-splitter transformation (J2J with |r| 2 = \t\ 2 = 1/2. Moreover, the field an 
after propagation is still described by a thermal statistics, since the Gaussian statistics and the factorization property 
P)l of the fourth-order correlation function are preserved by a linear unitary transformation. However, the coherence 
properties of the field an change upon propagation, as described by the well-known Van Cittert-Zernike theorem (see, 

e.g., ^3)- In particular, in the far-field plane where H (x, x') = —jj cxp ^— ^jyx ■ x'^j , they are described by a second 

order correlation function 

(cf H {x)a H {x')) otfdyf dy 'e>%(**-*'*')r (y, y ') . (22) 

By assuming translational spatial invariance of the source, the far-field correlation function Q22JI is proportional to 
S(x — x '), as can be easily verified by substituting Eq. © into Eq. (|22l) . However, when this unrealistic assumption is 
removed, the finite transverse size of the source ws has the effect that the correlation length of the function (|22|l is also 
finite and is inversely proportional to ws, as shown by the Van Cittert-Zernike theorem. Hence Z c ' oh oc represents 
the coherence length in the far-field plane. We can thus conclude that a high level of pixel-by-pixel correlation can 
be observed also in the far-field plane, provided that the size of the detection regions S is not too small with respect 
to i c ' oh , and the thermal beam is sufficiently intense (see the discussion in the Appendix). 

We remark that despite C can be made close to 1 by increasing the mean number of photons, it never reaches the 
quantum level, as shown by Eq. H20|) . 

For the entangled beams produced by PDC, spatial correlation is present both in the near and in the far field, with 
the ideal result (<5iV 2 } = 0, C = 1 in both planes [Tl| . In this case, the far field correlation is between symmetric 
pixels, and the coherence length in the far field is inversely proportional to the pump beam waist, which therefore in 



(21) 
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this context plays the same role as the source size w$ for the thermal beams. 

In we analysed the effect of replacing the pure PDC entangled state with two mixtures that exactly preserve the 
spatial signal-idler quantum correlations, either in the far or in the near field. It turned out that when considering 
the "far-field mixture" the pure state results could be exactly reproduced in the z = / configuration of Fig. [21 but 
no information about the image was present in the z = 2/ configuration. The converse is true considering the 
"near-field mixture" . This result is a consequence of the fact that the far-field intensity and the near-field intensity 
are non-commuting operators, and in the absence of quantum entanglement between the two beams they cannot 
simultaneously be correlated up to a perfect degree (this point is related to that raised in ,l6j). This led us to argue 
that only in the presence of quantum entanglement the whole set of results can be obtained by solely changing the 
setup in reference arm 2 0. 

However, nothing prevents two non-entangled beams to be correlated in both planes up to an imperfect degree. The 
two beams generated by splitting thermal light are actually imperfectly correlated both in the near and in the far-field; 
but by using intense thermal light, the classical intensity correlation is strong enough to reproduce qualitatively the 
results of both the z = f and the z = 2/ configuration. 



V. IMAGING PERFORMANCES IN THE CLASSICAL AND QUANTUM REGIMES 

A complete comparison of the performances in the classical and quantum regimes requires extended numerical 
investigations describing realistic thermal sources, which are outside the scope of this paper. 

However, some general remarks can be made concerning the key issue of the visibility of the information in the two 
regimes. The information about the object is retrieved by subtracting the background term (Ii(xi)) {I 2 {x2)) from 
the measured correlation function JfjJ, as indicated in Eq. Q. A measure of this visibility is given by evaluating the 
following quantity in relevant positions 

= G{x 1 ,x 2 ) = G{x u x 2 ) , 

{h{Si)h{x 2 )) (/i(fi))(/ 2 (f 2 )) + G(f 1 ,x 2 ) ' 1 dj 

with < V < 1. 

A first remark concerns the presence of (n(q))th in Eq. @ in place of U\{q)V2{— q) in Eq. As a consequence, 

in the thermal case G t h(xi, x 2 ) scales as {n(q))^ h . In the entangled case, GpdcK^i, ^2) scales as \U\{q)V2{— q)\ = 
(n(<f))pDc + ( n (<z))pDC' where (n(q))pDC = I V2 I 2 = |^i(<7)| 2 is the mean number of photons per mode in the PDC 
beams, and |£7i(<f)| 2 = 1 + | Vi(q)\ 2 (see, e.g., [Tl| ) - The difference between the two cases is immaterial when the mean 
photon number is large, while it emerges clearly in the small photon number regime ({n(q)) <C 1). Actually, in the 
thermal case the visibility does not exceed the value 1/2, whatever the value of (n(q))th, since Gth^i,^) scales in 
the same way as the background term. On the contrary, in the PDC case the visibility can approach the value 1 in the 
small photon number regime, since in this case the leading scale of Gpdc(^1j ^2) is {n(q))puc and this term becomes 
dominant with respect to the background (Ii(xi))(l2(x2)) °c ( n ( < ?)}pDC- Hence, in the regime of single photon pair 
detection the entangled case presents a much better visibility of the information with respect to classically correlated 
thermal beams (see also [l9jp. 

A second remark concerns the role of the temporal argument. Standard calculations show that the visibility scales 
as the ratio between the coherence time of the source T co h and the detection time (see also 0, EJ). This implies 
that conventional thermal sources, with very small coherence times, are not suitable for the schemes studied here. A 
suitable source should present a relatively long coherence time, as for example a sodium lamp, for which r co h ~ 10~ 10 s 
|15| . or the chaotic light produced by scattering a laser beam through a random medium (see, e.g., p^h 

As a special example of a thermal source, one can consider the signal field or the idler field generated by PDC. 
Fig. |3 shows the results of a numerical simulation j2(| for the reconstruction of the diffraction pattern of a double 
slit, in the scheme z = / of Fig- El It compares the use of the entangled signal and idler beams (curve (a)), and two 
classically correlated beams obtained by symmetrically splitting the signal beam (curve (b)). The parametric gain 
is such that (n(g))pDC ~ 750 at its maximum for (a), and (?i.(<7*))pdc ~ 1500 for (b), so that the the mean photon 
numbers of beams b±, 62 are approximately equal in the two simulations. From our simulations it clearly emerged 
that the number of pump shots necessary for reconstructing the diffraction pattern up to a desired accuracy is the 
same for both curves (a) and (b). Notice that Fig.[3]plots only G(xi,x 2 ), which contains the object information and 
represents the relevant part of the intensity correlation. This was obtained by subtracting a large background term 
(Jl(xi)) (12(^2)) from the correlation function of the intensities {Ii{x\)l2{x2)) ■ Consideration of these quantities (not 
shown in the figure) allowed us to calculate the visibility V, as defined by Eq.J2BJ, which turned out to be V ~ 0.05 
in both cases. Although this is a rather poor visibility, the crucial point is that the fringes shown in Fig. could be 
correctly reconstructed after a reasonable number of pump shots. 
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FIG. 3: Numerical simulation of the reconstruction of the diffraction pattern of a double slit in the scheme z = f of Fig. [5] 
G{xi,X2) is plotted versus X2 after 10 4 pump shots for (a) entangled signal/idler beams from PDC, (b) classically correlated 
beams by splitting the signal beam, (c) is the analytical result of Eq. 1141 . Parameters are those of a 4 mm /3-barium-borate 
crystal (l co h = 16.6 /im, r co h = 0.97 ps). The pump waist is 664 fim, and the pulse duration is 1.5 ps. The slits are 36 )im wide 
and slit separation is 122 /im. xq is defined as AqXf /(2tt). 




FIG. 4: Numerical simulation of the reconstruction of the image of a double slit in the scheme z = 2/ of Fig. [5] G(xi,X2) is 
plotted versus X2 after 10 4 shots for (a) entangled signal/idler beams from PDC, (b) classically correlated beams by splitting 
the idler beam, (c) is the analytical result of Eq. 1181 . Parameters as in Fig. [3] 



The simulations of Fig. [31 were repeated but changing the setup in the reference arm to the z = 2/ configuration. The 
results are shown in Fig. and confirm that the classical correlation can be used to reconstruct both the diffraction 
pattern and the image of the object by operating only in the optical setup in the reference arm 2. Also in this case 
the efficiency of the reconstruction is the same for curves (a) and (b). 

These examples clearly show that in the regime of high photon number the quantum and classical correlations offer 
similar performances. 

Another important aspect that emerges from these examples is that the classical and quantum imaging schemes 
apparently offer the same spatial resolution. This is most evident in the plot of Fig. ^ where the spatial resolution 
of both schemes is not good enough to reproduce the sharp details of the double slit image, but the reconstructed 
images are almost identical. As a matter of fact the analytical results - given by the general formulas ©-{TBI for 
the correlation function of intensity fluctuations, the results (|14fl - H15[) for the diffraction pattern reconstruction, and 
the results p& [l -p9 |) for the image reconstruction - show clearly that the spatial resolution of both schemes depends 
only on the spatial coherence properties of the sources. Provided that the spatial coherence properties of the classical 
source emulate those of the entangled source, there is no reason why the two schemes should offer different spatial 
resolutions. This requirement is not unrealistic at all, since the typical transverse coherence length of the entangled 
beams from PDC is on the order of tens of microns. An example of "thermal" light whose coherence properties can 
be engineered is offered by, e.g., chaotic radiation obtained by scattering laser light through random media |l8j. 

A final remark concerns the form of the analytical results of Eqs. (|14fl . (|15|) for the diffraction pattern reconstruction. 
Specifically, the diffraction pattern on the r.h.s. of Eq. <|15|l depends on the sum x\ + X2 of the positions of the pixels 
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in the detection planes of beam 1 and 2. This feature was exploited in the experiment of Ref. [2(| (see also [U), which 
was performed by registering coincidence counts of pairs of photons generated by PDC in a configuration similar to 
the z = f scheme of Fig. 27]. The interference fringes from a double slit were observed by scanning the pixel 
detectors in the two beams together, i.e. with x\ — x%. A halving of the period of the interference fringes with respect 
to those observed by illuminating the object with coherent light was then observed. This effect was claimed to be 
a consequence of the entanglement of the two-photon state. This effect is evident from Eq. Ijl5|l . where by setting 

X\ + 2^2 = x the PDC correlation function gives 



f ^2x|j^ . By inspecting Eq. I|14|) . we notice that the same effect 



could in principle be observed in the scheme that uses the classically correlated thermal beams, provided that the 
pixels in the detection planes of the two beams are scanned symmetrically, i.e. setting x\ — —x^ |28[ . Therefore, 
the discussion is open whether quantum entanglement has a crucial role in this observed halving of the period of 
the interference fringes, or whether this is just a consequence of the spatial correlation of the two beams and of the 
particular detection scheme used. 



VI. CONCLUSIONS 



In conclusion we have suggested a way of producing classically correlated beams suitable for correlated imaging. 
These beams are the outcome of mixing an intense thermal beam with the vacuum state on a beam splitter. We 
have shown a deep analogy between the use of entangled beams originating from PDC and these classically correlated 
beams. The analogy arises because of the similar structure of Eq. 10 and Eq. which implies that the outcomes 
of correlation measurements with the classical source can emulate those obtained with the entangled beams, for any 
choice of the test and reference arm optical setups. This holds provided the spatial coherence properties of the source 
are chosen to mimic the marginal statistics of the individual PDC beams. 

This analogy relies on the high level of spatial correlation that exists between the two beams emerging from the 
beam splitter as a consequence of the large excess noise of the thermal input. A key point is that a pixel by pixel 
correlation is present not only in the near-field plane immediately after the beam splitter but also in the far-field 
plane. Although the correlation is limited by shot noise, it allows to reconstruct both the image and the diffraction 
pattern of an object by only acting on the optical setup of the reference arm. This was hitherto thought of as a feature 
that truly required entanglement between the two beams 0. 

We have investigated the imaging performances of the quantum and classical schemes. Both our analytical results 
and a specific numerical example show that the spatial resolution limitations of the two schemes have similar origins, 
namely the finite transverse coherence length of the light. Thus, also in this respect the classical beams mimic the 
results of the entangled beams. 

On the other hand, in the small photon number regime a definite advantage of the quantum configuration is 
represented by a better visibility, as it was already recognized in other contexts (see, e.g., p^l. Thus, in imaging 
schemes where the visibility represents a crucial issue one should take this into account. However, as the number of 
photons per mode becomes large this advantage disappears and the visibility tends to be the same for the quantum 
and the classical source. This result suggests that the peculiar difference between the use of the two kind of sources is 
not given by the entanglement, but rather by the possibility of working in the photon counting regime in the quantum 
case. 

A further advantage of the PDC source may lie in the possibility of using a fraction of the source laser as a reference 
field in order to perform balanced homodyne detection. We will show in a future theoretical work that the homodyne 
scheme makes it possible to perform phase-sensitive correlated imaging with a high degree of visibility even in the large 
photon number regime [24[ • Another advantage of using the PDC source relies on the possibility of multi- wavelength 
imaging, as mentioned in the introduction. 

Our results implies that it is possible to perform coherent imaging without spatial coherence by using thermal light 
in combination with a beam splitter. This is reminiscent of the Hanbury-Brown and Twiss interferometric method for 
determining the stellar diameter [l3l | , as well as of detecting the fringes arising from interference of two independent 
thermal sources [23j. However, here we define a technique to achieve a fully coherent imaging of the object through 
correlation measurements with a high degree of flexibility. Moreover, since the required correlation is classical, a high 
quantum efficiency of the detectors is not necessary. 

The results of this paper suggest new experiments of correlated imaging with thermal light, which could provide a 
useful test bed for the imaging experiments with quantum entangled sources. Experiments like that simulated by the 
numerics of Sec. fVl may open new possibilities, offered by the combination of the correlated imaging from entangled 
beams and from classically correlated beams. 
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Appendix 

This Appendix investigates how the degree of spatial correlation between the beams obtained by splitting thermal 
light depends on the size of the pixels used to detect the light, and on the mean number of thermal photons. 

As shown by Eq. I|21|) a high level of spatial correlation can be present between the beams b\ , 62 after the beam 
splitter as a consequence of of a high level of excess noise in the thermal beam a. The relevant quantity to consider 
is the ratio (: SNf :)/ (N\), where (: SNf :) is defined by (5Nf) = (Ni) + (: SNf :), and represents the noise in excess 
with respect to the coherent state level. 

By using the beam splitter transformation and the factorization property (J3J of the fourth-order thermal 
correlation function, we get 

(: 5Nf :) = |r| 4 f dx f dx' \T(x- x')f . (24) 
Js Js 

where T(x — x') is the thermal correlation function defined by Eq. J2J). On the other hand, the mean number of 
photons detected over the pixel is 



(Ai) = |r| 2 Sr(0) cx |r| 2 -^n max , (25) 

'coh 

where we introduced the parameter 

"max = (n(q = 0)} th (26) 

that represents the mean number of photon in the most intense mode q = in the spectrum. Moreover, we used 
r(0) = J j0yz{n(q))th oc l/^ 2 oh n max , where the proportionality constant depends on the actual shape of the spectrum. 
When the pixel size is much smaller than ^ co h, the correlation function in Eq. (|24|l is approximately constant over the 
integration regions 5, so that 

(: 5Nl :) -» |r| 4 S 2 r(0) 2 = (A : ) 2 . (27) 

Hence, in this limit of a small detection region the noise takes the form of the usual single-mode result for thermal 
light (SNl) = (JVi) + (Ai) 2 . However, as shown by Eq. (gSJ the mean number of photons is also small in this limit, 
so that the excess noise and the correlation tend to be both small. 

On the other side, when the detection regions grow much larger than the coherence area the excess noise does not 
scale any more with the square of the mean number of detected photons. In the limit of large detection regions the 
r.h.s. of Eq. (|2"4l can be approximated as 



trv .plane 



2 



dq 



(: 6N( :) - |rf S / d£ r(0 = |r| 4 S / —L-(n(q)f th 



« |r| 4 ^-T4««|r| a fw<^i>. (28) 

Voh 

For a large detection region the excess-noise scales proportionally to the mean number of detected photons, so that 
the ratio (: SNf :)/{Nx) reaches a limiting value, which depends on n m ax- 

We have investigated more in detail the role of the size of the detection region in the degree of correlation by assuming 
square pixels with sides of size A, as well as a Gaussian correlation function r(£) = r(0)e~^! 2 / i <=° h , which is a good 
approximation of a smooth correlation function decaying on a distance l co h- By substituting this in Eq. I|24() and 
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using the result of Eq. 1(21(1 , we get the two relevant asymptotic behaviours of the degree of correlation 1(21(1 : 



C 
C 



1 



S 2 n n 



4tt 1 + ±6 2 n r 



41 



for S < 1 



for S — > oo, 



(29) 
(30) 



where 6 — A/Z co h is the ratio of the pixel size to the coherence length. The curves in Figure [3] show the general 
behaviour of the degree of correlation C as a function of the pixel size for different values of the mean photon number 
per mode n max . The consequence of this figure is that when the input thermal beam is intense enough, a high degree of 
spatial correlation can be achieved. This is provided that the pixel size is not too small with respect to the coherence 
length describing the decay of the thermal spatial correlation function. 




FIG. 5: Degree of spatial correlation C between two identical detection regions of the beams obtained by splitting thermal 
light, as a function of the ratio 5 between the pixel size and the coherence length. n max is the mean photon number in the most 
intense mode. C = 1 represents the maximum degree of correlation. 



Another relevant parameter that characterizes the spatial correlation between two beams is the visibility of the 
correlation, which in the spirit of Eq. 121(1 we define here as 

_ WjVg) 

Vs - ImW (31) 

This definition is analogous to that of Ea. 1(23(1 . where, however, a small pixel was implicitly considered so that no 
integration over the pixel area was performed. This quantity is easily calculated in the same manner as above. Figure^] 
plots the visibility of the spatial correlation between the thermal beams as a function of the pixel size scaled to the 
coherence length. This plot was obtained under the same assumption of Fig. The visibility turns out not to depend 
on the mean thermal photon number n max , so that only one curve is plotted. 




FIG. 6: Visibility Vs of the spatial correlation between two identical detection regions of the beams obtained by splitting 
thermal light, as a function of the ratio S between the pixel size and the coherence length. 



Unlike the degree of correlation, the best visibility is obtained when the detection pixel is not large with respect to 
the coherence length. The main conclusion of this Appendix is therefore: in order to achieve both a good visibility as 
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well as a high degree of spatial correlation, the best choice is a detection pixel with a size approximately equal to the 
coherence area. 
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There is unfortunately a misprint in Eq. (8) of Ref. Q, so that the argument of the function T is X2 instead of — X2- 
In our numerical simulations we dropped the translational invariance in both space and time of the input beam, and thus 
we also took into account the temporal variable (which was ignored in the analytical treatment). The simulations were done 
in 2D+1 dimensions, i.e. including one transverse dimension (along the walk-off direction) as well as the time dimension, 
and propagating along the ^-direction of the crystal. The calculated intensities were integrated over time since the pulse 
length was much shorter than the response time of any available detectors with spatial resolution. See [Tl| for more details. 
[27] The scheme of the experiments of Ref. I2H is slightly different from the z = / scheme of Fig. |2] and corresponds to 
having an object in both arms. However, this only has the effect that the Fourier transform of the square of the object 
transmission function is observed instead of the Fourier transform of the object transmission function, which is irrelevant 
when the object is a double slit 

[28] One can easily imagine a setup where in arm 1 there is an additional optical element that flips the beam transverse 
distribution, so that x\ can be scanned together with $2 as in the experiments of [20^1 and pl|. 



